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ABSTRACT 

We study the time evolution of the wave function of a particle bound by an attractive 8- 
function potential when it is subjected to time dependent variations of the binding strength 
(parametric excitation). The simplicity of this model permits certain nonperturbative 
calculations to be carried out analytically both in one and three dimensions. Thus the 
survival probability of bound state |6>(t)| 2 , following a pulse of strength r and duration t, 
behaves as |6>(t)| 2 — |#(oo)| 2 ~ t _a , with both 9(oo) and a depending on r. On the other 
hand a sequence of strong short pulses produces an exponential decay over an intermediate 
time scale. 

* Research Supported by AFOSR Grant F49620-98-1-0207 
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1. Introduction 

While there has been much progress in our understanding of the processes leading to 
the ionization of atoms and/or the dissociation of molecules subjected to time dependent 
fields, the mathematical difficulties presented are such that there are no explicitely solvable 
models for transitions form a bound into the continuum [1-15]. This motivates us to 
investigate here the ionization probability of a particle bound by an attractive point 5- 
function potential in one dimension [5,6,15] and a spherically smeared out 5-f unction in 
three dimensions. We obtain explicit expressions for the ionization probability and for 
the energy distribution of the ejected electrons for certain time dependent parametric 
excitations, i.e. when we suddenly change the value of the coupling constant for a time 
interval t. For such changes the survival probability of the bound state shows no regime 
of exponential decay but approaches its asymptotic value as a power law [11-13]. The 
situation is different for periodic forcing with short pulses which is also treated here and 
more generally in [16]. The survival probabilities now include intermediate exponential 
regimes followed by power law asymptotics. 



General formulation 



We consider first the one-dimensional system with unperturbed Hamiltonian [5,6,15] 



Hn = 



2m dx 2 



gS(x), g>0, 



-oo < X < oo. 



(1) 



H has a single bound state 

u b (p,x) = v^e- p|x| , P=^9 (2) 

with energy —E = —fkvo = —ti 2 p 2 /2m and a continuous uniform spectrum on the positive 
real line, with generalized eigenfunctions 

u(p, k, x) = -i= ( e lkx ^—-S kx \ J , -oo < fc < oo, (3) 

and energies h 2 k 2 /2m (with multiplicity two for ^ 0). Here u b is normalized to 1 and 
u(k, x) to 5(k — k'). 

Beginning at some initial time, say t = 0, a perturbing potential V(x, t) = —R(t)S(x) 
is applied to the system, i.e. we change the parameter g in H , 

g^g + R{t), t>0. (4) 

We note here that the matrix elements, |(itb|y|/c)| 2 = R 2 (t)-^k 2 /(p 2 + k 2 ) which vanishes 
as k — > and approaches R 2 (t)-^ as \k\ — > oo. This implies in particular that the integral 
of the transition matrix over all k is infinite. 

To solve the time dependent Schrodinger equation, 

ih^f^ = H iP(x 1 t)-R(t)5(x)iP(x 1 t) 1 t>0 (5) 

we expand ij){x,t) for t > in the complete set of functions u: 

o roc 2 

ip(x,t) = 6(t)u b (p,x)e i ^ t + e(k,t)u(p,k 1 x)e- i ^ t dk 7 t>0 (6) 

J — oo 



and monitor the evolution of 6{t) and Q(k,t) starting from the initial bound state 

9(0) = 1, 0(fc,O) = o. 

The ionization probability at time t caused by a pulse, which coincides with R(t') for 
t' < t and vanishes for t' > t, is given by 

/oo 
|e(M)| 2 ^, (7) 
-oo 

while |6>(£)| 2 is the survival probability. 

This model can be extended to a three dimensional shell-like delta function potential. 
The Hamiltonian, 

H Q = - — A -gS(r - a), a > 0, r = |r|, r e M 3 (8) 

ill 

has bound states with angular momentum / for all I = 0, 1, ... such that / < mga/h 2 — 1/2. 
The time dependent perturbation is now of the form V(r, t) = —R(t)5(r — a). The results 
for three dimensions, which are similar to those in one dimension are described in Section 
5. This follows calculations of 6{t) and Q(k,t) in one dimension. 

2. Integral equation for the 1-d case 

Using the orthonormality of the eigenfunctions (2), (3) and substituting (6) into (5) 
yields the following set of equations for the time dependent amplitudes at t > 0, 



where 



T(t) = 



1 



i\k\ 



V2n J-oo P + i\k 



■ e -^ {p2+k2 ^@{k,t)dk 



R(t) (10) 



determines both 9 and 6: 



n Jo 



Q(k,t) = 



\k\ 



f T{t'y^ p2+k2 »'dt'. 

Jo 



(11) 
(12) 



'2tt{p - i\k\)h , 

Substituting (11) and (12) into (10) yields an integral equation for T(t) which using 
dimensionless variables obtained by setting h = 2m = g/2 = 1 (implying p = 1, uiq = 1) 
yields 

6{t) = l + 2i f Y{t')dt', S(k,t) = , / Y(t')e l{1+k2)t 'dt\ (13) 

Jo V 7T 1 - i\k\ Jo 

where Y(t) is to be found from the integral equation 



Y(t) = rj(t) |l + jf [2i + M(* - t')]y(t')cft' 



and r7(t) = R(t)/g. The function M(s) in (14), is given by 

2i f°° u 2 e - JS (i+ w 2 ) 



M( S ) = - / — - 
^ Jo 1 



M(s) behaves as 



M(s) 



-j=^e- %s + 0(s" 5 / 2 ), whens^oo, 
-4±L -i + O^ 1 / 2 ), whens^O. 

V27TS V ' 



(14) 



(15) 



(16) 



3. Ionization by a rectangular pulse 



We study perturbations having the form of a step function, rj(t) = r for t > 0. The 
calculation of P at any time t will then correspond to the ionization probability caused 
by a pulse of amplitude r and duration t. Substituting the rj(x) into (14) and taking the 
Laplace transform we find, c.f. [16] 



r 



Y(s) = Y{x)e- S *dx = — (17) 

s — r(i + yis — 1) 



o 



where 

%Vis~^T>0. (18) 
Using the inverse transform Eq.(13) has the form 



r fioo+a e st _ ^ 



6{t) = 1 + - / -: — p T ^ ds, a>0. (19) 

7T J_ ioo+a s[s - r(i + VlS - 1)\ 



To evaluate (19) we make a cut in the complex plane of s along the imaginary axis 
from — too to —i. In the left half plane bounded by the left side of the cut and the vertical 
line from — % + a to +zoo + a the integrand in (19) is analytic except for a simple pole at 
s = ir(r + 2) when r + 1 > 0. There are no poles if r < —1, i.e. when the coefficient of the 
5-function is positive and the potential for t > represents repulsion. The integral along 
the left half-circle of infinitely large radius is clearly zero, therefore one may rewrite (19) 



as 



9{t) = - / - ds + 2 

tv J ss- (i + Vis^l)} 



r + 1 + \r + II 



s[s - (i + Vis~^T)] " {r + 2f 
6 



ir(r+2)t 

& 1 



with counterclockwise integration around the cut. Straightforward manipulations with the 
integral term allow us to write 6 finally in the form 

e(t) - ir! r e ~ l{1+u2)t u Hu + 2 r+1 + ' r + V ^ (20) 

The integral in (20) can be expressed in terms of Fresnel's functions and the dependence of 
the ionization probability on r is shown in Figs.l and 2 where it seen that it is monotone 
for r < — 1 but not for r > 1 so we can have "atomic stabilization" [14,15]. 

Using (12) and (17) one can calculate \<S)(k,t)\ 2 which gives for t > r the energy 
distribution of electrons kicked out from the bound state by a pulse of duration r. We find 
in the original units 



'2p i\k\(p — q) J (q + \q\)e l ^ T e~ l ^^ T 



7T (p-i\k\) I (q 2 + k 2 )(p + q) (p + i\k\)(q-i 



p + q 



^ T Erfc ( i±l.f^i\ _ i?i e ^ r Erfc ( l+l./Vr 

p 2 + k 2 \ 2 V m I p 2 + k 2 [ 2 V m 



+ 



(p 2 + k 2 ){q 2 + k 2 ) y 2 V m J { 

where q = (1 + r)p and Erfc(z) denotes the probability integral. For large k, |G(/c,r)| 2 

decays like k~ 4 giving a very long tail to the energy distribution of the emitted electrons. 

In Fig. 3 we plot \B(k, r)| 2 vs. k for several values of r when r = —1, i.e. when the pulse 

just destroys the attractive interaction. It is seen that the longer the pulse the more peaked 

is the distribution with the maximum moving towards small values of k. 
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The total energy of the electrons ejected by the pulse is given by 



fc2 roo 



r)\ 2 dk, t>r. 



For measurements made outside the range of the potential this energy will be the same as 
the kinetic energy of the emitted electrons. An analytical evaluation of E(t) yields a very 
long and not very illuminating formula. Instead we present in Fig. 4 a numerical plot of 
E(t) for r = — 1. When the pulse length r — > oo, E(t) approaches the value 



E(oo) = E 



r + 1 + 1 



l + 2|r + l| + 2(r + l + |r + l|) 



(r + l)(r + 3) 
(r + 2) 2 



which increases linearly with |r| when |r| — > oo: 



£(oo) -> 2£ |2r + |r||. 



Attractive long pulses, r > 0, thus give three times as much energy to the ejected electrons 
than do the repulsive ones, r < 0. This is shown in Fig. 5 where it is seen that E(oo) is 



monotone for both positive and negative r. 



We note that a rectangular pulse perturbation is a special case of a sudden jump from 
the initial Hamiltonian Ho to a new time-independent Hamiltonian Hi which in return 
jumps to H when the perturbation ends. The amplitudes 9, and 6 can thus also be 
calculated by projecting (twice) a new state onto the old one which is just the evaluation 



of overlap integrals. The Laplace method, which gives 9{t) for general tine dependence of 



R(t) in the form 

i rioo+a st -| 

= 1 + - / -y(s)ds, (22) 

shows that |6>(t)| 2 will have an exponential decay for t — > oo only when F(s) has poles in 
the left half of the complex plane s. 



Power law decay 

When the pulse length t goes to infinity the integral in (20) vanishes and in the limit 
we have 

| 9(oo)P = /l6(r+l)V(r + 2)', ifr>-l, 

[0, if r < -1 v ' 

It is seen from (23) that any two very long pulses (at least one of them must be repulsive) 
produce the same ionization if their amplitudes r and r' satisfy the relation - + ^ = — 1. 

For large t the asymptotics of the integral term in (20) can be easily found. Using 
contour integration we can rewrite the integral as 

i r°° y 2 e~y 2 

(r + ijwi Jo FwFw+?i y 

and integrate by expanding the integrand in powers of y 2 /t. Let us first study the case 
r = — 1, which corresponds to the perturbation removing the potential and making the 
electron evolve for t > like a free particle. The decay of the bound state in this case is 
rather slow: 

|^)| 2 = l + 0(t- 2 ), r=-l. (24) 



When both t and t\r + 1| are large we get 

= 2 r + l + |r+ll <tr(r+2) ^— = e - rt + 0(t" 5 / 2 ). (25) 

For the survival probability of the bound state we have 



1 W(oo)\ 2 + 8r 2 cos[(r+l) 2 t] . { _ 1 
l^)| 2 ^<| 4 (r+1)(r+2)2tV ^' (26) 

if r < -1. 



(r+l) 4 7rt 3 ' 

Thus for r < — 1, when the evolution takes place with a repulsive 5-function, the approach 
to zero of |6>(t)| 2 is like t~ 3 , compared to the t~ x decay given in (24) for the free evolution, 
see Fig. 2. Note that the coefficient of t~ 3 becomes independent of r for \r\ » 1. For 
r > — 1 the approach of |6>(t)| 2 to its nonvanishing asymptotic value is oscillatory with an 
envelope which decays like t -3 / 2 . These oscillations are very rapid for large r (Fig.l), but 
their amplitude is small, of order 1/r. These asymptotic power law decays are in agreement 
with general results for the decay of initially localized states, c.f. [5-8]. 

4. Ionization by periodic short pulses 

The behavior of P(t) for short pulses of duration t « 1, is very different for cases 
when a = ry/i is large or small comparing with 1. Writing P(t) = P(t,a) we analyze 
Eq.(20) and have in the case of a single pulse 

P(t,a)=4 A /^( 2a2 / 3 fora«l, (2?) 
V ' V 7T \ 1 for a » 1. V 7 

We turn now to the survival probability when we bombard our system with a whole 

train of short pulses of duration r << 1 repeated periodically with period a ~ 1. Using 
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(13) yields 

n 

0(na + t) = 6 n = 1 + 2i J k, ®n = 1 if n < 0, 
where we have defined 



(28) 



fc=0 



J™ = / (na + r- x) m Y(x)dx. 

J no 



(29) 



By integrating Eq.(13) for 9{t) in t from t = no to t = no + r and using (16) we obtain 



m+l 



r _1 <C = r 

m + l 



n-l 



l + ^(2z + M n _ fc )4° 

A;=0 

_J jm+1 _| 3km 

m + l n V^(2m + 3) 



+ \l-k m J™ +1 / 2 + 



(30) 



J™+ 3 / 2 +0(r m + 2 ), 



where M n = M(no), k m = ^/nr(m + 1)/T(m + §). The inequality \6 n \ < 1 implies 
| J? | < 1 and therefore by integration by parts we get 



I Tin ^ 

\ J k I < T ■ 



(31) 



Let us eliminate in (30) the term J™ +1 / 2 by using (30) with m — > m + 1/2 which gives 

% T m+3/2' 



m+l 



r _1 <c = ! 7 , 

'm + l 



+ 2/c m r 



m + l 



+ 



7r 2m + 3 
r/c m /c m _|_i/2 



n-l 



7T 



l + ^(2z + M n _ fc )J fc ° 

k=0 

J™ +1 +0(r m+2 ). 



+ 



(32) 



Combining (31) with the estimate |M(s)| < a/7t/2s 3 we have an upper bound on the sum 
in (32) in the form | X]fc=o M n -kJ® \ < 2.4cr -3 / 2 max je [ 0)n _i] | J°|. Treating the amplitude 
r as a quantity of order of unity one can immediately improve the upper bounds (31) for 



J™ to: 



\' J n ' 



(33) 
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Eq.(33) allows us to rewrite Eq.(32) for m = as a simple reccurence 



n-l 



Jl = p l + 2z^4° +r 2 /n(r), 



(34) 



fc=0 



where for a > 1, r < 1 we have |/ n (T)| < 7 uniformly in n and 



p = rr 



l 4rV?(l + i) 



3V27T 



Starting with Jq = p ~ r one can find successively J° using (34). The terms of such a 
sequence will be close to the corresponding terms of the solution of the simplified equations 



n-l 



J n = p 1 + 2i Jk , 



(35) 



fc=0 



till \ J n \ » t 2 and | J n \ » \p\nr 2 . It is easy to solve (35) to find 



J n = p(l + 2ip) n , 



(36) 



and therefore \J n \ |p|e n7 where 7 = 8r 2 r 3/ ' 2 /3v / 27r << 1. Our condition for 
approximating J° by J n has now the form 



e~ ni » nr 2 . 



(37) 



Using J n and (28) we find 



Q(t) » exp ( " 7 + 2zrT t ] . 



(38) 



if the duration of train of pulses t = no is not too long and satisfies (37). One can obtain 

from (37) that the decay of survival probality \9(x) \ 2 up to a value p is accurately described 
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by (38) if yfr < 2r 2 ^/Ji/ ln/j -1 ; for fx = 0.01, r « 1 this gives yfr < 0.04 and a train of 
about 300 pulses. For shorter r the train can be longer and the ionization more complete. 

For longer trains of perturbation the term X]fc=o M n -kJ® in (30) cannot be ignored 
and it makes the eventual decay slower with strong oscillations due to interference with the 
eigenfrequency. We see that in the exponentially decaying regime the survival probality 
is independent of a. This is very different from the case where the time dependent 
r](t) = rsmuit considered in [16]. In that case the exponential decay depends strongly 
on to. In our case r — * which means that 77(f) will contain all ranges of frequencies. 

5. Three dimensional model 

The Hamiltonian (8) has eigenfunctions in the continuum spectrum 

*i, m (k,r) = Y ljm (6,<p)R l (k,r), r > 0, < 9 < tt, < <p < 2tt, 
where the radial functions are 

Ri(k,r) = A l ^J l+1/2 (kr) + 

fk f 0, if r < a, 

A l r \ f QJ ;+1/2 (te)< / #a)^ ) 1/ #r) - H^ 1/2 (kr)H^ 1/2 (ka)] r > a. ^> 

The dimensionless parameters A\ normalize Ri(r) to a 5-function, 

Mk) = {1 + QnJ l+1/2 (ka)N l+1/2 (ka) + Q 2 n 2 Jf +1/2 {ka)[jf +1/2 {ka) + W, 2 +1/2 (fca)]}~ 1/2 , 

and the notations for normalized spherical harmonics Y^ m and Bessel functions are the 

usual ones. The energy corresponding ipi,m(k,r) is h 2 k 2 /2m. 
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The parameter Q = 2mga/h 2 plays a crucial role for the existence of the bound states, 

QK l+1/2 (p l a)I l+1/2 (p l a) = 1 (41) 

is the equation for the energy —h 2 p 2 /2m of all bound /-states (they are of different axial 
symmetry). The left side of (41) is a monotonically decreasing function of its argument 
7 = pia and it is equal to Q(2l + l) -1 when 7 = 0, therefore 

Q > 21 + 1 

is the condition to have the bound states for all V <l. 

The radial normalized eigenfunctions can be written in the form 

R\{r) = ^{Y 1/2{ r ] <K ( MK < \ if C" a ' ( 42 ) 
y/f \Ii + i/2{Pia)Ki +1/2 {pir)/K l+1/2 {p l a) 1 r > a, 

where 

y/2K l+1 / 2 (pia) 

l = — 7 ============================= 

y/1 - PiaKi +1/2 (pia)[Ii_ 1/2 (pia) + Ii +3/2 (pia)} 

and /, K are the modified Bessel functions. 

There are no transitions between states of different angular symmetry if both the 
potential in (8) and perturbation V(t, r) are central. For simplicity we consider our three- 
dimensional model with Q > 1 in in the s-state. Dropping the index / = equation (41) 
for the energy of the bound state —ti 2 p 2 /2m is 

2ap 



Q = 



^ g — 2ap ' 
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The eigenfunctions (42) of the bound and the continuum states are respectively 



*&(r) 



P 



1/2 



f sinhpr, 



if r < a, 



V7r(e 2pa - 1 - 2pa) \ e p(r a) sinhpa, r > a, 
2 -1 / 2 fsin/cr, if r < a, 



7rrWl - Q 



sin 2fca _i_ f)2 sin 2 



r v> sin^ ka \ « m At - Q 8 -^ sin fc(r - a), r > a. 

ka T fc 2 a 2 



(43) 
(44) 



Assuming that the particle is in the bound state ^(r) at t = and the perturbation 
has the form V(r,t) = —R(t)g5(r — a), we use the method of projections which was 
described in Section 4 to find the ionization probability induced by the rectangular pulses 
R(t) = rg for t > 0. After the end of pulse at t = r we have for 6>(r) an equation similar 
to (21), 



9{r) = 



4pq 



( e 2 P a _ 1 _ 2pa)(e 2 i a - 1 - 2qa) 



e (p+q)a p e (p-q)a _ q e (q-p)a 



p + q 



(45) 



8p 



[(pa — Qi) sinhpa + pa cosh pa]' 



2 r°° 



na* 



L 



e 2m sin rea 



(p 2 + P) 2 (i-Qi^ + Q 2 ^) 



a7c, 



where q is the solution of Eq.(41) with Qi = (1 + r)Q instead of Q, (g gives the energy 
of new bound state). If Qi < 1 the first term in (45) vanishes, otherwise the square of its 
absolute value represents the probability 1 — P(oo) of the electron to remain in the bound 
state when r — > oo. Using the dimensionless time uj t — > t the asymptotics of the decaying 
term in (45) when t — > oo is 



0(t) = 0(oo) - v 7 ^ 



[(pa — Qi) sinhpa + pa cosh pa] 2 1 + i 



+ 0(t" 5/2 ), t^oo, (46) 
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or 

\ft(i\\ 2 « Iflfn^l 2 -U i 

0(t" 3 ), if Qi < 1. 

The dimensionality as one can see changes the character of asymptotics only of the free 
evolution (t~ 3 / 2 vs t -1 / 2 ). An interesting case is Qi = 1, when the perturbed Hamiltonian 
has a "zero energy bound state". The asymptotic behavior of 6{t) is now given by 

m = 4[(po - gi) sinhpq + paooshpa] 2 (l - Q + Q(t - 3/2)> ^ ^ (4?) 

p 2 a 2 V27rt 

which has the same character as for the free decay in the one-dimensional model. 

In three dimensions the same technique as that used in Section 3 allows us to derive 
a one dimensional integral equations similar to (14) for each pair of quantum numbers 

/, m < I: 



T hm {t) = R(t)a 2 



itf(a)0j, m (O) + -J K l (t- t')T Um {t')dt' 



(48) 



which determines the evolution. In particular, the amplitude of the bound state develops 
in time as 

^(t)=^(0) + i^ J\ m {t)ti! . (49) 
The function Ki in Eq.(48), 

POO 

K0) = [R^(a)] 2 + / |fl,(fc,a)|V <R ( fc3+p ?W 2m dfc, (50) 
Jo 

is independent of the quantum number m. Each spherical harmonic evolves autonomously 

and if 6>^ m was zero at t = it does not change for our perturbation. Though the kernel 
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of Eq.(48) even for / = 0, 



4psinh 2 pa 2 f 00 e lh 13 sin 2 ka 

° ( } = a\e^ - 1 - 2pa) + ^ ^ j _ gsin2fca + g 2 sh^fca ^ 51 ) 



cannot be expressed in terms of standard functions numerical calculations are quite feasible. 



6. Concluding remarks 

Some general features of our results with possible implications for realistic systems 
include: 

a) The ionization probability approaches its asymptotic value as t~ 3 / 2 if the electron 
can be bound in the perturbed state, goes to zero as t~ x if the perturbation makes the 
electron a free particle, and as t~ 3 when the perturbation converts the attractive well into 
a repulsive one. 

b) A finite train of periodically repeated short pulses makes the survival probability of 
the bound state decay exponentially without oscillations. When the frequency of repetition 
is comparable with the eigenfrequency of the bound state or is lower the decay scales in 
such a way that only the total number of pulses is important. 

c) The three dimensional potential gives a similar behavior of the ionization. The free 
evolution in one dimension corresponds here to a marginal situation with the "zero-energy" 
bound state. 
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FIGURE CAPTIONS 



Fig. 1. The survival probability of bound state following the imposition of rectangular 
pulses of duration t for different relative amplitudes r > —1, see Eqs.(23),(24). 

Fig. 2. The normalized survival probability following the imposition of repulsive 
rectangular pulses, r < —1. Plots of \6{t) /9 as {t)\ 2 vs. t, where |6> as (t)| 2 is 4/nt for r = — 1 
and (r/r + l) 4 /nt 3 if r < -1, see Eqs.(24),(26). 

Fig3. Plot of |0(&;, t)| 2 which represents the energy distribution of electrons kicked 
out by the rectangular pulses of duration t, see Eq.(21). 

Fig 4. The total energy of states of electrons in the continuum spectrum ejected by 
the rectangular pulse of duration t for r = — 1. 

Fig. 5. Plot of the electron kinetic energy vs. the amplitude of very long rectangular 
pulses. 
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